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Abstract 

The number of two-dimensional percolation clusters whose external hulls 
enclose an area greater than A, in a system of area fi, behaves at the critical 
point as C^l/A for large A, where C = l/87r\/3- Here we show that away 
from the critical point this factor is multiplied by a scaling function that is 
asymptotically proportional to a simple exponential exp{—A/A*) where A* 
scales as \p—Pc\~'^'^- The fit is better than for Kunz and Souillard sub-critical 
scaling, which would predict the asymptotic behavior exp{—{A/A*)'^^^) 
where D = 91/48 is the fractal dimension. 



1 Percolation 

In the percolation model, one considers a disordered system, constructed 
(for example) by taking a regular lattice and diluting it by making only a 
fraction p of the sites or bonds conducting (occupied). One is concerned 
with the appearance of long-range connectivity in the system - the perco- 
lation transition - and the behavior near the point of that transition. For 
definiteness, we consider mainly bond percolation on a square lattice, and 
define 

s = the number of wetted sites of a cluster. (1) 
Isolated sites with no bonds attached correspond to s = 1. 
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Bond percolation on a square lattice, showing the appearance of a 
percolating path from the top of the system to the bottom 

Following are pictures of a single system where p is slowly increased 
(bonds are added) and clusters are merged (as in the Newman Ziff algo- 
rithm) . 
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2 Traditional percolation distribution 

What is the size distribution of clusters in percolation? According to scaling 
theory (e.g., StaufFer & Aharony 1994), the number of clusters of size s per 
site of the lattice behaves as 

ns{p) = cis-^f{c2{p-Pc)sn (2) 

where r and a, are universal exponents, and f{x) is a universal function. 
However ci and C2 are non-universal, but vary with the lattice and percola- 
tion type being considered. Here /(O) = 1 and f{x) — > as x — > ±00. The 
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scaling relation can also be written as 

ns{p) = c,s-^F±{s/s*) (3) 

where s* = is a typical size and scales as \p — Pc\~^^'^ ■ If a ratio of 
moments such as M3/M2 is used to define s*, then s* will have a coefficient 
that differs above and below pc, and that difference will have to be included 
in the definition of -^±(2;) in terms of f{z). In any case, F±{x) is universal, 
but the non- universal coeffiecient ci remains. 

Indeed, can never directly be put into a completely universal form 
because it concerns the size s which is a lattice-level quantity and by its 
very nature non-universal. 



3 Universal form of the size distribution 

We have shown (Ziff, Lorenz, & Kleban 1999) that, at the critical point Pc, 
the distribution of cluster size can be written in a completely universal form: 

A'"(clusters whose enclosed area > ^) ~ — ^ (4) 

in a system of total area CI, where C is a universal constant which for 2d 
percolation is given by (Cardy Ziff, 2003) 

C = -^-= « 0.0229720373 (5) 

Results have also been found for the more general Potts model. 
This result can also be phrased: 

• The number of clusters per unit area A whose enclosed area is greater 
than A equals a universal number C. 



4 Proof of (5) 

Prom the definition of Ug we have 

Hsds ~ ns^-'' (6) 

Each cluster has one external hull which encloses an area 

A ~ r2 ~ (7) 
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where D is the fractal dimension of the cluster and r is a length scale. Thus, 
s ~ and 

iV(clusters whose enclosed area 

= n/A (8) 



by virtue of the hyperscaling relation 

i = r-l (9) 

Furthermore, the constant C must be universal, since A^(area > A) is only a 
property of the larger clusters, or, in other words, a property of the universal 
percolation fractal. 



5 Derivation of the exact value of C 

Eq. @ implies that the average depth, defined as the number of outer hulls 
crossed in going from a point to the edge of the system, scales as 

Clog(OMo) (10) 

where Aq is a lower cutoff of the area (~ the lattice spacing). Transforming 
conformally from an annulus to a rectangle using z — > logz, we find that 
this corresponds to a rectangular system with SvrC clusters crossing per unit 
length. But Cardy has shown that that number is l/\/3, which implies 



as given above. 



6 Numerical demonstration 

We generate "rooted" hulls by carrying out a hull generating walk on a 
lattice, and stopping when the cluster closes or when an upper cutoff is 
reached. Using a (virtual) lattice of 65536 x 65536, and a cutoff of 2^^^ = 
4194304 steps, we add a weight to each walk of a factor of l/(no. steps), and 
thus find an unbiased estimate of A^(area = A) as long as A <cutoff. We 
bin logarithmically, 

N{A, 2A) = (hulls w/enclosed area G {A, 2A))/VL (12) 
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N{A,2A) = C 



which should behave asymptoticahy for large A as 

1 _ 1S\ - _^ 
^~ 2A) ~ 2A 

or 2AN{A,2A) = C (or 2C for both interior and exterior hulls): 

2A N(hull area in range (A,2A)) vs. Log_2 A 

0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 



Levels out to the asymptotic value 
of 1/(4 pi sqrt(3)) = 0.0459... 



(13) 
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7 Off-critical behavior 

Away from pc, we expect 

CO 

iV(hulls whose enclosed area > ^) ~ —G{A/A*) (14) 

where A* = typical hull area and G{x) is a scaling function that should be 
universal. One expects 

A* ^ ^—^ (15) 

where ^ is the correlation length and v = 4/3 in 2d, since the area of a 
percolating cluster is not a fractal (even though its boundary is a fractal). 

Here we plot the same quantity as before but at p = 0.495. The quantity 
AN {A, 2A) for both internal and external hulls approaches the value of 2C 
for log A ~ 5 but then deviates for larger A. 
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A N(A,2A) vs. log_2 A 
p = 0.495 
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8 Average Area 

We can define an average area by averaging over all clusters (normalizing 

by ny. 



so that 



(A) = ■i^A2Ar(area = A) 

- / A'^N{area =A)dA 
" Jo 
1 /"^ 

= - / 2AiV(area > A) dA 

/>oo 

= 2C G{A/A*)dA 
Jo 

/•oo 

= 2CA* / G{x)dx (16) 

JO 

(A) cx ^* (17) 
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Note that we do not use ^ A Ar(area= A) to define the average area, because 
that quantity scales as biA*/Ao; where Aq is a lower cutoff (area of basic 
cell of the lattice). In fact, that quantity gives the average depth of crossed 
loops. 

9 Numerical measurement of average area 

We have verified that (A) indeed scales as \p — Pc\~'^^ ■ 
p {A) logarithmic slope 



0.4 39.786 

0.42 70.351 -2.554 

0.44 148.060 -2.587 

0.455 313.931 -2.612 

0.46 427.332 -2.618 

0.465 607.705 -2.620 

0.47 912.729 -2.639 

0.475 1478.133 -2.628 

0.48 2670.927 -2.646 

0.485 5730.272 -2.652 

0.49 16841.970 -2.659 



Table. Average area size, per lattice site (both internal and external 
hulls) as a function of p. The logarithmic slope A log A/ A log (]9c — p), 
clearly approaches a value consistent with —8/3. 

Here it is plotted on a log-log scale: 



14 



Mean area <A> vs. (pc - p) on a log-log scale 




10 Higher-order behavior of the mean area 

Assuming an exponent of 8/3, we can find the higher-order behavior by 
plotting (pc — p)^^^{A) vs. {pc — pY , and varying x until we find a straight 
line, which occurs for x = 4/3: 
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(Pc - P)^^^ <A> VS (Pc - P) 
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which imphes 



(Pc - P)' 



{A) = aQ{pc-p) ^''^ + ai{pc - p) 



-4/3 



(18) 



11 Asymptotic behavior of G{x] 



For large x, we find G{x) decays as e 
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Exponential decay of scaling function 
for internal plus external hulls 
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y = -0.12980X + 0.43911 
= 0.99996 



A/<A> 



Because the precise definition of A* is arbitrary up to a constant, we can 
fix that constant by requiring that G{x) is a simple exponential for large x: 



which implies that 
so that 



G{x) ~ e X > 1 

A* = (^)/0.130 

Jo 



2C 



(x)dx = 0.130 



(19) 
(20) 

(21) 



12 Comparison to sub-critical percolation mass 

For s 3> s*, the distribution of mass s in subcritical percolation is believed 
to behave as (Kunz and Souillard 1978): 

ns{p) - e-'^'/'* (22) 
s 

• This is evidently consistent with our result only if s ~ A (instead of 
s ~ A^^"^ = ^91/96^ [qj. subcritical clusters 
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• Note also that Ug apparently has a more complicated form than A^, 
since changes its leading power term from s~'^ for s ^ s*, to s^^ 
for s ^ s*, while N always has the same leading behavior of for 
all A. 



13 Amplitude ratio of average areas 

According to the universal of G{x), the amplitude ratio of 



(23) 



is also universal. (Note - there is only one A* in the system, used for scaling 
for both internal and external clusters). The measurements suggest a value 
of about 180 for this ratio (note higher statistical error as p ^ Pc)'- 



Amplitude ratio of <A^'^>/<A'"*> 
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14 Alternate form of scaling 

We can also write as 

CO 

iV(area > A) = ^g{±{A/A*)^/'^''] 



(24) 
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where l/2v = 3/8 and it refers to external (+) and internal hulls (— ). 
We write the argument of g in this way because then it is proportional to 
{p ~ Pc)^^'^'^ and therefore g{z) should be a single analytic function of z for 
both positive and negative as in Q for the size distribution. Here are 
our results for p = 0.495 for the scaling function g{x): 

g(x) 



1 *^ 

3^ 
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(A/<A>)'^{3/8) 



15 Conclusions 

We have shown 
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iV(hulls whose enclosed area > ^) ~ ^G{A/A*) (25) 

with G{x) a unique universal function satisfying 

G{x) = 1 X = (26) 

G{x) ~ ce"^ X > 1 (27) 
where c ^ exp(0.439) = 1.55 and A* = (^)/0.130. 



16 Appendix: Scaling for N(area = A) 

For A^^ = number of hulls whose area equals A, we can write 
Na = ^F^{A/A*) = ^-f{c,{v-Pc)A"^l 



(28) 
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where l/2v = 3/8 and f{x) = F{x^'^). This is not a universal form of the 
distribution both because it involves clusters of exactly area A, and also 
because the non-universal constant C3 shows up in the argument of / (in the 
second form). It is interesting however as the counterpoint of the scaling 
form of the percolation mass that is also analytic at p = 0. 

Here we keep internal and external hulls separate. By duality, we have 
the symmetry, 

t'''\x) =f''\-x) (29) 

which allows us to the entire curve of /(x) by making measurements at 
p < Pc only, and using the internal hulls for —x and the external hulls for 
X. Here are the results for p = 0.44, where we plot A'^Na vs. (p — Pc)A^/^), 
analogous to ((11)) : 
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17 Note 

An inconsistent factor of 4 was used in the definition of the area, and this 
error may not have been consistently corrected in all of the plots. 
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